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PREFACE 

This  report  was  prepared  by  Dr.  Yin-Chao  Yen,  Research  Physical  Scientist,  of  the  Geotechnical 
Research  Branch,  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions  Research  and  Engi¬ 
neering  Laboratory.  The  work  was  funded  under  DA  Project  4A161 101 A91D,  In-House  Labora¬ 
tory  Independent  Research,  Work  Unit  350,  Growth  of  Phosphorus  Smoke  in  Winter. 

Technical  review  of  the  report  was  performed  by  Dr.  Yoshisuke  Nakano  and  Dr.  Virgil  Lunar- 
dini  of  CRREL. 
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At  Surface  area  of  droplet 
cd  Heat  capacity 

D  Diffusion  coefficient  of  water  vapor  in  air 
d  Diameter 

/  Molar  flux  of  vapor 

/a  Molar  flux  of  air 

K  Defined  as  pdj/pd 

k  Thermal  conductivity  of  air 

L  Latent  heat  of  condensation 

Molecular  weight  of  water 
m  Mass  of  droplet 

n  Molar  concentration 

Pd  Vapor  pressure  over  droplet 

Pt  Saturation  vapor  pressure,  saturation  humidity 
PiQ  Equilibrium  vapor  pressure  over  plane  water  surface 
4h  Percent  relative  humidity 
^rhs  Percent  of  PRH  at  droplet  surface 
r0,r,  Defined  in  eq  25 
f , ,  r  i  Time  derivative  of  r,  and  r2 ,  respectively 

r, ,  r2  Laplace  transform  of  r,  and  r2 ,  respectively 
r  Spatial  coordinate 

R  Radius,  gas  law  constant 

Sd  Ratio  of  PJPto 

t  Time 

rc  Characteristic  time  defined  as  R\jD 
T  Temperature 

u  Transformation  variable  defined  as  u  -  rpw, 

u  Integration  factor,  defined  in  eq  2 1 

Vt  Molar  volume  of  water 

V2  Molar  volume  of  air 

x,  Mole  fraction  of  water 

Subscripts 

a  Air 

«>  Bulk 

d  Droplet 

ds  Droplet  surface 

vwz  Water  vapor  in  air 
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RHS  Relative  humidity  at  droplet  surface 

s  Saturation,  surface 

0  Initial  (also  superscript) 

v  Vapor 


Greek  fetters 

a  Thermal  diffusivity 

0  Defined  in  eq  21 

p  Density  of  vapor  or  droplet 

p'  Density  of  solution  droplet 

co  Defined  asp® a /pds 

r  Defined  as  t/tc 

6  Absolute  temperature 
o'  Surface  tension  of  solution  droplet 
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AEROSOL  GROWTH  IN  A 
COLD  ENVIRONMENT 


Yin-Chao  Yen 


GENERAL  BACKGROUND 

Aerosols  are  commonly  formed  either  by  the  conversion  of  gases  to  particulate  matter  or  by 
the  disintegration  of  liquids  or  solids.  Aerosols  may  also  result  from  the  resuspension  of  powdered 
material  or  the  break-up  of  agglomerates.  Formation  from  the  gas  phase  tends  to  produce  much 
finer  particles  than  disintegration  processes.  Dust,  smoke,  fumes,  haze  and  mist  are  all  terms  in 
common  use,  with  somewhat  different  popular  meanings.  Dust  usually  refers  to  solid  particles 
produced  by  disintegration  processes,  while  smoke  and  fume  particles  are  generally  smaller  and 
formed  from  the  gas  phase.  Mists  are  composed  of  liquid  droplets.  Generally,  the  generic  term 
aerosol  is  used  to  describe  all  such  systems  of  small  particles  suspended  in  air  or  some  other  gas¬ 
eous  medium. 

Atmospheric  aerosols  are  usually  composed  of  a  mixture  of  particulate  materials  originating 
from  various  sources,  including  gas-to-particle  conversion  processes.  The  sources  are  a)  the  com¬ 
bustion  of  fuel  oil  in  power  plants  and  refineries,  including  sulfur  dioxide  and  fuel  oil  flyash,  and 
b)  automobile  exhaust,  including  particulate  matter  composed  of  lead  bromochloride  and  tarry 
organic  matter.  The  gas  phase  includes  nitrogen  oxides,  organic  vapors,  and  sulfur  dioxide.  A 
small  percentage  of  these  gases  is  converted  to  nitrates  resulting  from  the  oxidation  of  NO  and 
N02 ,  oxygenated  organic  compounds,  and  sulfates.  This  conversion  process  increases  the  total 
aerosol  concentration  and  contributes  to  degradation  of  visibility  significantly. 

Soil  dust  that  becomes  airborne  as  a  result  of  the  wind  and  of  surface  activities  associated  with 
transportation,  construction,  and  agriculture  also  contributes  to  the  aerosol  concentration.  But 
because  of  its  relatively  larger  size  (>  1  ftm)  it  tends  to  settle  out,  although  it  is  continuously  re¬ 
introduced  by  surface  disturbance. 

Aerosols  may  affect  human  health,  visibility  and  climate.  Particle  size,  concentration,  and 
chemical  composition  are  usually  the  most  significant  factors  determining  such  effects.  Light 
scattering  by  ..mall  particles  is  also  a  sensitive  function  of  size.  The  maximum  scattering  efficiency 
corresponds  to  particle  size  of  the  same  order  as  the  wavelength  of  the  incident  light.  Processes 
leading  to  the  accumulation  of  particles  in  that  size  range  produce  the  most  severe  visibility  degra¬ 
dation.  Such  processes  depend  on  the  chemical  composition  of  the  system,  which  also  determines 
the  refractive  index,  another  factor  in  li^it  scattering  and  absorption.  The  nucleating  efficiency 
of  small  particles,  i.e.  the  ability  to  form  clouds,  is  expressed  in  terms  of  the  rate  of  condensation 
per  unit  mass  of  aerosol.  Particles  smaller  than  the  critical  size  determined  by  the  Kelvin  effect 
(because  of  increased  vapor  pressure  over  a  small  particle)  do  not  grow.  Large  particles  have  a 
relatively  small  amount  of  surface  area  per  unit  mass,  which  accounts  for  their  reduced  nucleating 
efficiency. 


A  complete  treatise  of  the  fundamentals  of  aerosol  behavior,  i.e.  aerosol  characterization,  trans¬ 
port  properties,  deposition  by  convective  diffusion,  internal  deposition,  light  scattering  and  vis¬ 
ibility,  experimental  techniques,  processes  of  collision  and  coagulation,  gas-to-particle  conversion, 
and  thermodynamic  properties,  can  be  found  in  the  book  by  Friedlander  ( 1 977).  This  book  is 
aimed  at  enhancing  the  understanding  of  the  science  of  particulate  pollution  and  applying  it  to 
practical  problems,  such  as:  How  are  aerosols  formed  at  pollution  sources?  How  can  particles  be 
removed  from  gaseous  emissions  to  prevent  them  from  becoming  an  air  pollution  problem?  And 
how  can  air  quality  be  related  to  emission  sources  to  devise  effective  pollution  control  strategies? 

A  comprehensive  review  on  cloud  droplet  growth  was  presented  by  Carstens  and  Carter  (1974). 
Briefly,  the  process  of  growth  involves  the  transport  of  mass  (vapor)  toward  the  drop,  release  of 
latent  heat  at  the  drop  surface,  subsequent  heating  of  the  drop,  and,  as  a  consequence,  thermal 
energy  transport  away  from  it.  The  continuity  of  water  vapor  concentration,  n  (mole/volume), 
outside  the  drop  is 

'■'-F  0) 

where  /  is  the  molar  flux  of  vapor,  given  by 

-  (n  +  /!.)  D  r 

/=  (1  (7*  +  a'  (x)  Vlnfl+y—-  7«  (2) 

where  jr  =  mole  fraction  of  vapor 

n  =  molar  concentration  of  vapor 
D  =  diffusion  coefficient  of  water  vapor  in  air 
T  =  temperature  field 

/a  and  nt  =  molar  flux  and  concentration  of  air 
a’(x)  =  thermal  diffusion  factor. 

Under  normal  atmospheric  conditions,  further  simplifications  can  be  made  and  the  flux  and  the 
corresponding  heat  conduction  become 

05‘p- 57  (3) 


«V7-=if 

with  three  of  the  boundary  conditions  as 
P  (°°.  t)  =  p„  (t) 

7’(«,0=7’-(0 


4tr/?*  \kj£-  +  LD  1^1  =cd 

l  dr  9rJ  r=R  d  dt 


where  p  *  vapor  density 

D  and  o  =  diffusion  coefficient  of  water  vapor  in  air  and  thermal  diffusivity  of  air -water 
vapor  mixture,  or  simply  that  of  air 


p„  and  r„  =  bulk  values  of  vapor  density  and  temperature  (which  can  be  functions  of  time) 
R  =  drop  radius 
k  =  thermal  conductivity  of  gas 
L  =  latent  heat  of  condensation 
cd  =  heat  capacity  of  drop. 

With  the  assumption  that  the  transients  appropriate  to  the  diffusion  problems  are  small  enough 
so  that  the  steady-state  profiles  “follow”  the  outward  motion  of  the  drop  surface,  the  flux  equa¬ 
tions  reduce  to 

V2p  =  0  (8) 


V2r=0.  (9) 

I 

Boundary  conditions  5  and  6  remain  the  same,  but  condition  7  reduces  to 

k  I  +LD  ^  I  =0.  (10) 

dr  lr=R  dr  lr=R 

Based  on  the  material  balance  on  the  drop  surface  the  outward  deep  growth  rate  is 

—  =  —  *1  (11) 

dt  pd  dr  I  r=  R  uu 

where  pd  is  the  liquid  density. 

In  general,  two  somewhat  different  growth  treatments,  both  based  on  the  above  assumptions, 
exist  in  the  meteorological  literature.  The  most  common  approach  assumes  that  the  transport  is 
entirely  controlled  by  the  diffusion  of  mass  and  conduction  of  heat,  while  the  alternate  approach 
posits  the  possibility  of  additional  control  exerted  at  the  liquid/vapor  interface.  In  what  is  per¬ 
haps  the  most  common  treatment  of  drop  growth,  it  is  assumed  that,  to  a  good  approximation, 
thermal  equilibrium  between  vapor  density  and  temperature  exists  at  the  drop  surface.  Further¬ 
more,  the  dependence  of  thermal  equilibrium  on  both  drop  curvature  and  dissolved  salts  must  also 
be  taken  into  account.  The  alternative  approach  is  that  the  drop  growth  is  controlled  by  diffusion 
and  surface  kinetics.  The  vapor  molecules  do  not  always  stick  to  the  surface  upon  striking  (or  are 
inhibited  from  evaporation),  and  consideration  of  coefficients  of  condensation,  evaporation  and 
thermal  accommodation  is  necessary.  The  thermal  equilibrium  conditions  used  in  diffusion-con¬ 
trolled  growth  are  replaced  by  mass  and  heat  fluxes.  This  is  commonly  accomplished  by  equating 
the  molecular  and  energy  fluxes  (at  the  drop  surface)  determined  from  uniform  gas  kinetics  to 
those  calculated  directly  from  Fick’s  and  Fourier’s  laws. 

Most  of  the  works  reviewed  and  summarized  by  Carstens  and  Carter  (1974)  are  concerned  with 
theories  of  drop  growth  in  naturally  occurring  aerosols  aimed  at  meteorological  applications.  Little 
work  has  been  reported  on  man-made  chemical  or  nonchemical  aerosols.  These  aerosols  are  rather 
short-lived  in  nature  and  can  serve  as  a  means  for  camouflage  or  as  an  obstacle  to  target  recogni¬ 
tion. 

Rubel  (1978)  seems  to  be  the  only  investigator  to  have  conducted  a  detailed  analysis  of  phos¬ 
phorus  smoke.  In  his  work,  the  equilibrium  vapor  pressure  over  the  surface  of  a  droplet  containing 
an  aqueous  solution  of  volatile  solvent  and  nonvolatile  solute  was  uniquely  determined  by  Kelvin 
and  solute  effects.  Expressions  for  the  saturation  ratio  over  the  surface  of  a  phosphoric  acid  drop¬ 
let  as  a  function  of  droplet  diameter  and  moles  of  acid  are  presented  as 


P,  4  o'  Mx 
In  = _ -In 


R  Tdp 


r. *  — l 

L  nd3-6n2  {V2  +  F,)J 
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4.55+ - - — —  ,0. 
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577  >  x,  >  0 


(13) 


where 


Pd  =  vapor  pressure  over  droplet 

PiQ  =  equilibrium  vapor  pressure  over  plane  surface  of  pure  water 
o'  and  p'  =  surface  tension  and  density  of  solution  droplet 
A/j  =  molecular  weight  of  water 
R  =  universal  gas  constant 
T  =  temperature  of  system 


r/,  and  V2 
x. 


moles  of  acid 

molar  volume  of  water  and  acid 
mole  fraction  of  water  in  acid  solution. 


Equations  12  and  13  relate  the  saturation  ratio,  i.e.  Sd  =  PJPio,  over  a  phosphoric  acid  drop¬ 
let  to  its  fundamental  properties,  its  diameter,  and  the  number  of  moles  of  acid  («2)  contained  in 
it.  In  both  expressions,  the  first  term  on  the  right  represents  a  vapor  pressure  increase  due  to  the 
Kelvin  effect  and  the  second  term  represents  a  vapor  pressure  reduction  due  to  the  solute  effect. 
Rubel  (1978)  reported  further  that  for  droplets  containing  dilute  solutions,  eq  12  can  be  simpli¬ 
fied  as 


In 
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12  n2  K, 


R  Tp'd 


ltd3 


(14) 


Though  this  equation  is  not  applicable  to  highly  concentrated  solution  droplets,  it  nevertheless 
clearly  reveals  the  Kelvin  effect  (first  term  on  the  right)  and  the  solute  effect  (second  term  on  the 
right)  as  functions  of  droplet  diameter.  For  small  diameters,  the  solute  effect  dominates  the  Kel¬ 
vin  effect,  resulting  in  the  saturation  ratio  Sd  =  PJPso  being  less  than  1 .  As  the  diameter  increases 
and  the  solution  becomes  dilute,  the  two  effects  become  equal  (i.e.  Sd  =  1),  and  thereafter  the 
Kelvin  effect  dominates,  resulting  in  Sd  >  1 .  As  the  droplet  diameter  d  increases  to  infinity, 
both  effects  become  negligible,  and  the  system  thus  physically  represents  a  plane  surface  of  pure 
liquid  with  5d  =  1 . 

To  establish  a  relationship  between  the  relative  humidity  RH  and  the  droplet  diameter,  Rubel 
(1978)  further  assumed  that  diffusional  equilibrium  exists  between  the  surface  of  the  droplet  and 
the  environment,  and  obtained  predictive  relationships  between  RH  and  droplet  diameter  by 
replacing  Sd  =  PJPto  by  RH  in  eq  12  and  13.  He  computed  the  droplet  diameter  as  a  function 
of  RH  and  found  that  the  relative  increase  in  droplet  diameter  is  independent  of  the  initial  con¬ 
densation  nucleus  diameter  dc,  and  that  the  ratio  dKH/dc  (where  </RH  is  the  droplet  diameter  at 
that  specific  relative  humidity)  is  only  dependent  on  RH.  The  relative  increase  in  this  ratio  takes 
on  values  of  1 .5, 2.0  and  3.0  for  RH' s  of  70, 90  and  98%,  respectively,  irrespective  of  condensa¬ 
tion  nucleus  variations  from  0.30  to  2.40  pm.  The  results  clearly  demonstrated  that  the  solute 
effect  dominates  the  Kelvin  effect,  and  since  the  ratio  dRH/dc  remains  nearly  constant  for  specific 
RH  values,  all  droplets  achieve  the  same  degree  of  dilution  and  experience  the  same  relative  in¬ 
crease  in  diameter,  irrespective  of  initial  condensation  nucleus  size.  Assuming  diffusional  growth 
equilibrium  between  the  environment  and  the  droplet  surface,  Rubel  (1978)  concluded  that  the 
solution  effect  governs  the  properties  of  a  phosphoric  acid  droplet  in  phosphorus  smoke,  and  that 
the  greatest  proportion  of  growth  occurs  in  the  region  defined  by  relative  humidity  greater  than  80%. 
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Rubel  also  considered  this  problem  from  the  kinetics  point  of  view  in  order  to  verify  the  de¬ 
velopment  of  the  relationship  between  relative  humidity  and  droplet  diameter,  i.e.  the  adoption 
of  diffusional  equilibrium  between  the  environment  and  the  H3P04  droplet.  He  reported  that  the 
relaxation  times  (characteristic  times  for  approach  to  equilibrium)  are  extremely  low  for  relative 
humidity  from  10  to  98%— in  the  range  from  10-3  to  0.90  second.  In  other  words,  for  the  case 
of  an  unsaturated  environment,  acid  drops  attain  equilibrium  with  the  environment  almost  instan¬ 
taneously,  and  the  assertion  of  diffusional  equilibrium  under  tactical  conditions  is  proved  to  be 
valid. 

Rubel’s  phosphorus  smoke  work  can  be  represented  by  the  following  sequence  of  reactions: 
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lM 


P4  +j  Cl  -*  P40m 

4(»)  ^(*)  4  10(V)(») 


P4°‘0(v)(J)  +  6H2O(v)->4H3PO4(8) 


(V) 


’(solution) 
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Cragin  (1982)  conducted  some  snow  sample  chemical  analyses  during  CRREL’s  SNOW-ONE 
field  experiment  at  the  Camp  Ethan  Allen  Training  Center,  Burlington,  Vermont,  in  connection 
with  smoke  tests  to  evaluate  the  effectiveness  of  instrumentation  used  to  measure  transmissivity 
and  other  optical  properties.  The  concentration  of  the  smoke  material  deposited  on  the  snow 
downwind  can  be  expressed,  in  general,  by  an  exponential  decay  function,  i.e.  C- A  exp  (-Bx). 

The  constants  .4  and  B  are  found  to  be  dependent  on  the  nature  of  the  chemical  smoke  as  well  as 
the  wind  velocity.  Since  the  scale  of  the  test  was  very  limited,  the  results  reported  by  Cragin  can 
not  be  considered  conclusive;  they  can  only  indicate  the  trend  of  the  smoke  particle  disposition 
process. 

Since  chemical  smoke  growth  rate  has  not  been  exclusively  studied  under  a  variety  of  condi¬ 
tions,  this  research  is  aimed  at  deriving  an  expression  approximately  representing  smoke  aerosol 
growth  as  a  function  of  time  and  humidity  at  subfreezing  temperatures.  To  accomplish  this,  assump¬ 
tions  will  be  made  to  circumvent  the  difficulty  arising  from  the  continuous  dilution  of  the  conden¬ 
sation  nuclei  by  the  process  of  chemisorption. 


THE  PROBLEM 

The  analysis  presented  here  concerns  the  rate  of  growth  following  the  instantaneous  formation 
of  condensation  nuclei.  A  mass  balance  at  the  drop  surface  results  in 


—  =DA 
dt  UA* 


dp, 


Wt 

dr  1  r=R 


05) 


where  m  -  mass  of  droplet 

D  =  water  vapor  diffusion  coefficient 
A ,  =  spherical  drop  surface 
Pw»  =  water  vapor  density  in  air 
t  and  r  =  respectively  time  and  spatial  coordinate. 

Equation  1 5  can  be  rewritten  in  terms  of  droplet  radius  as 


a  dR_n  *Pw. 

fi*dT-Dir 


r=R 


(16) 
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with  an  initial  condition  of 


R  =  R0  at  t  =  0  . 


Considering  spherical  symmetry,  the  diffusion  of  water  vapor  through  the  air  mass  to  the  enlarging 
aerosol  can  be  written  as 


3pw»  „  I 

^2Pwa  2 

3Pwa\ 

(17) 

St 

^  dr2  r 

— j 

with  boundary  conditions  of 

Pwa  “  Pds 

at  r  =  R 

P  =  p® 

^wa  ''wi 

at  t-  0 

(18) 

where  pds  is  the  water  vapor  density  at  the  drop  surface  and  p®  a  is  the  water  vapor  density  at  time 
zero. 

With  the  substitution  of  u  =  rpwa,  eq  17  and  18  can  be  replaced  as 


du 

dr 


=  D 


d2u 

dr2 


(19) 


and 


u  -  pdJ/?  at  r=R 

u  =  rp®  m  at  t  -  0  .  (20) 


METHOD  OF  SOLUTION 

Multiplying  both  sides  by  the  factor  exp  (-(Jr)  and  integrating  with  respect  to  r  over  the  range 
from  R  to  and  with  the  introduction  of 

u  =  f°°  u  exp  (-0r)  dr  (21) 

Jr 

eq  19  after  some  simplification  becomes 

exp  (D  fPt)  ~  [exp  (-D02t)u)  =  -  (p  ~  +  DQu  +  u  |  =r  exp  (-0R).  (22) 

Multiplying  eq  22  by  the  factor  exp  (- DfPt ),  integrating  from  0  to  <*>,  and  using  the  definition  of 
u,u  =  rpwt  and  eq  1 6,  eq  22  can  be  transformed  into 

P~f  0  +^o)«p(-^o)  =  /“  ^pd,+0pd,0*  +  pdK  f  +  *Pd,  f) 


exp  (-0 R)  exp  (- D02t )  dt . 


(23) 


For  large,  pure  water  drops,  the  value  of  pdj  will  be  equivalent  to  the  saturation  vapor  density  over 
a  plane  water  surface.  However,  for  small  water  droplets  containing  dissolved  salts,  the  value  of 
pdj  will  be  dependent  on  the  curvature  of  the  droplet  and  the  solute  concentration.  To  alleviate 
these  complications,  a  study  of  the  effect  of  pdt  on  the  aerosol  growth  rate  was  carried  out,  with 
the  assumption  of  constant  relative  humidity  maintained  over  the  droplet  surface. 

Defining  the  dimensionless  parameters  K  and  w  as 


and  noting  the  fact  that  the  value  of  K  (for  the  case  of  water)  is  in  the  range  of  10-8  to  10-7,  and 
with  the  introduction  of  s  =  D  ,  eq  23  becomes 

“*[f+  V7*°j,x4  i! *<]-/■  ♦>/*«♦**  f] 

exp  (-ccR)  exp  (-$f)  dt .  (24) 

A  series  solution  in  terms  of  the  small  dimensionless  parameter  K  can  be  expressed  as 

R  =  r„  +  Krx  +  K2r2  +  K3r3  +  K\  +  ....  (25) 

The  detailed  calculations  in  the  expression  of  rn  (n  -  0, 1 , 2, ...)  are  given  in  the  Appendix.  De¬ 
fining  a  characteristic  time  fc  *  Rq/D  and  a  dimensionless  time  t  =  t/tc,  the  expression  for  aerosol 
radius  up  to  K 4  as  functions  of  oj  and  r  can  be  represented  by 

-ij  -K(cj-  1 )  (1 .1284  r*4  +  t)  +  K3  (w-1)2  (0.7184  t*  + 

+  0.8634  t  -  0.3716  W  -  0.5  i2)  +  K3  (p-  l)3  (0.5555  t 

-  0.3451  r  +1.2628  t3!2  -  0.6265  r2  -  1.2412  t*/2  +  0.5  t3) 

+  K*  ((w-  l)3  (-0.8165  r1/2  -  3.2035  r  +  0.2394  r3/2 

+  0.0131  r2  +  0.0376  t5/2  -  0  1667  r3)  +  (w-  1)4  (-0.5507  t1/2 

-  3.9272  t  +  5.2030  r3/2  -  1.6403  r2  -  0.5600  t*/2  +  0.8910  r3 

-  1.1660  t7/2  -  0.5  r4)]  .  (26) 

RESULTS  AND  DISCUSSION 


To  evaluate  droplet  growth  rate  under  subfreezing  conditions,  environmental  temperatures  of 
0°,  -5°  and  -10°C  were  considered.  At  each  temperature,  four  levels  of  water  vapor  saturation 
(i.e.  relative  humidity  PRH  *  25,  50, 75  and  100% Pt)  were  used.  To  account  for  the  effect  of 


v  X. 


Table  I .  Values  of  K  and  to 


1.1448  2.2895  3.4343  4.5790 

0.9290  0.4579  0.6869  0.9198 
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0.4579 

0.9158 

1.3749 

1.8316 

4.856 

9.684 

14.580 

19.400 

0.6869 

1.3737 

2.0606 

2.7474 

7.264 

14.530 

21.790 

29.050 
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2.7474 

3.6632 
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29.050 

38.740 

^RH  (mm  H8) 

0.7908 

1.5815 

2.3723 

3.1630 

P dug  -  20%  Pou 

0.1582 

0.3163 

0.4745 

0.6326 

0.4745 

0.9492 

5.112 

10.230 

PRHS  “  60%  ^RH 

0.6326 

1.2656 

1.8978 

2.5304 

1CX108 

6.814 

1 3.630 
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Id 
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rsrio  i>n 
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^RHS  =  P RH 
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0.3224 
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dissolved  salt  in  lowering  the  vapor  density  on  the  droplet  surface,  once  again  four  levels  of  vapor 
saturation  at  the  droplet  surface  at  each/*RH  (i.e.  />RHS  =  20, 40, 60  and  80%  ofFRH)  were  in¬ 
vestigated.  For  a  given  value  of  the  corresponding  values  of  K  and  a>  were  computed  from 
K  =  pdj/pd  and  w  =  ,/pds  respectively  and  are  shown  in  Table  1 .  Once  the  values  of  K  and  w 
are  known  the  quantity  (R/R0  -  1 )  can  be  computed  from  eq  26. 

Figures  1-3  show  the  effect  of  the  environmental  air  temperature  and  relative  humidity  on  the 
aerosol  growth  rate.  It  is  evident  that  the  growth  rate  increases  with  a  decrease  in  /’rhs  and  with 
an  increase  in  t  for  t  up  to  a  value  of  106 .  It  is  also  clearly  demonstrated  that  the  growth  rate  in¬ 
creases  with  the  relative  humidity.  Since  the  log  (R/R0  -  1 )  vs  log  (r)  plots  are  approximately 
linear  with  r  for  r  greater  than  1 0,  cross  plots  of  (R/R0  -  1 )  vs  PRHS  with  PRH  as  parameter  are 
shown  in  Figure  4  for  Tt  ~  0°,  -5°  and  -1 0°C  and  at  a  r  value  of  1 04 .  For  each  temperature,  an 
identical  relationship  can  be  formulated  and  expressed  as 

(*T  )■””(£)  000 A"s)  (27> 

where  a25  is  the  slope  of  {R/R0  -  1)  vsPRHS  at  PRH  =  25%  Pt,  and  on  is  at  any  PRH  greater  than 
25% Pj.  Table  2  summarizes  the  on  values  at  three  levels  of  temperature  and  at  PRH  equal  to  25, 
50, 75  and  100%  Ps. 

From  these  tabulated  values  of  on,  it  can  be  seen  that  the  ratio  of  on /o2 s  can  be  expressed  in 
terms  of  PRH  (percent  of  Ps)/2S  for  all  the  temperatures.  A  plot  of  on  vs  6  (temperature  in  Kel¬ 
vins)  using  PRH  as  parameter  is  shown  in  Figure  5.  It  can  be  seen  that  the  dependence  of  on  on 
temperature  is  identical  for  all  the  values  of  PRH  covered  in  this  study.  The  variation  of  an  with 
temperature  can  be  represented  by 


a„  =  2 .2883  x  10-14  exp  (0.07368  0 ) . 


Replacing  o2s  by  eq  28  and  the  ratio  of  on/o2S  by  PRH/25  and  plotting 


ft-)  A 


,2883x  10-14  exp  (0.07368  0)  (100-PRHS) 


(28) 


vs  t,  all  the  computed  values  of  (R/R0  -  1 )  shown  in  Figures  1  -3  can  be  brought  to  a  single  line, 
and  the  values  of  (R/R0  -  1)  can  be  well  represented  by 

(iT"])  =  2491 7*  10-18  exp(O.O737  0)  (100-PRHS)  (r)0  9890.  (29) 

Once  the  size  of  the  condensation  nuclei,  i.e.  R0,  the  environmental  temperature,  the  value 
PRH,  and  the  assumed  value  of  PRHS  are  given,  eq  29  can  be  used  to  compute  the  aerosol  growth 
rate.  It  can  be  seen  that  for  the  same  r,  the  radius  of  the  aerosol  increases  with  temperature  ex¬ 
ponentially  and  is  directly  proportional  to  the  factors  of  PRH  and  (100  -PRHS). 

From  the  definition  of  r  =  Dt/R J,  for  an  initial  value  of  R0  =  0.15  pm,  a  diffusion  coefficient 
of  water  vapor  through  air  D-  0.26  cm2/s,  0  =  273K ,PRH  =  50%Ps,and  an  average  value  of 
PR hs  =  50%  PRH,  it  takes  r  =  0.01  s  to  obtain  a  value  of  R/R0  -  1  =  0.3343.  However,  based  on 
calculations  given  by  Rubel  (1978)  for  a  hydroscopic  phosphoric  acid  droplet  of  initial  nucleus 
radius  0.1 5  pm,  it  took  only  about  10"4  s  to  reach  an  equilibrium  size  of  R-  0.195  pm  in  a  50% 
relative  humidity  environment.  It  is  about  100  times  faster  than  the  case  of  a  pure  water  droplet 
(or  ice).  The  discrepancy  may  be  attributed  to  the  fact  that  the  rate  of  aerosol  growth  depends 
on  the  availability  of  water  moisture  in  the  air  (i.e.  absolute  humidity).  In  the  present  analysis, 
a  subfreezing  environment  is  considered,  and  thus  there  is  a  much  lower  absolute  humidity  than 
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Prhs<*Pr«1 


PRMS  ^PRH* 


a.  Ta  =  0°C,  Ps  =  4.579  mm  Hg. 


b.  Ta  =  -5°C,  P,  -  3.163  mm  Hg. 


* 


Ta  -  -10°C,  P,  =  2.149  mm  Hg 


Figure  4.  Variation  of  fR/R0  -  1 )  with  P^  HS  a 
t=  104  and  PRH  at  25,  50,  75  and  100%P,. 
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Figure  5.  Variation  of  on  with  6  at  PRH  *  25, 50,  Figure  6.  Variation  of  fR/R0  - 1 )/ 2.2883*  10~4  exp  (0.0737 
75  and  100 %  P4.  0)x  P RH/25  (100  -  ?RHSJ  with  r  for  all  ?RH,  PRHS  and 

temperature  values. 


Table  2.  Values  of  oa  as  functions  of  temperature 


*RH-*P, 


2S 

0.121X1 0"4 

0.086X1 0-4 

0.060X10 

SO 

0.242 

0.172 

0.120 

75 

0.363 

0.258 

0.180 

100 

0.484 

0.344 

0.240 

in  the  case  of  the  experimental  work  conducted  by  Rubel.  Furthermore,  the  processes  of  chemi¬ 
sorption  of  water  vapor  by  the  hydroscopic  nucleus  and  subsequently  the  dissolution  of  the  acid 
nucleus  may  have  greatly  enhanced  the  processes  of  transporting  water  vapor  to  the  surface  of  the 
enlarging  aerosols. 


CONCLUSIONS 


An  equation  for  transporting  water  vapor  through  an  air  mass  to  an  enlarging  water  droplet  has 
been  established  and  solved  with  the  method  of  Laplace  transformation.  The  series  solution  is  ex¬ 
pressed  in  terms  of  a  small  dimensionless  parameter  K  (as  shown  in  eq  25).  Droplet  radius  varia¬ 
tion  with  K  up  to  A4  has  been  computed.  The  results  covered  three  temperatures  and  four  levels 
of  relative  humidity  (25, 50, 75  and  100%  of  the  saturation  value).  A  simple  expression  relating  the 
instantaneous  droplet  radius  to  its  original  radius,  environmental  temperature,  relative  humidity, 
percent  of  relative  humidity  on  the  droplet  surface,  and  a  dimensionless  time  r  was  developed. 
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APPENDIX.  EVALUATION  OF  rn  s  IN  EQUATION  25 . 

Equation  25  is 

R0  =  r0  +  Krx  +  K2r2  +  K3r3  +  K4r4  +  ...  (25) 

with  initial  conditions  of 

ro  |t=o=/*o>rn  1»2*3  •••)•  (A1) 

Substituting  eq  25  into  eq  24  with  r0-  R0,  expanding  the  exponential  term  exp  (-s/s(D  R)  into 
a  series,  and  retaining  the  terms  up  to  K4,  we  have 

r(  / d(R0  +  Kr.  +  K2r2  +  K3r3  +  K4r 4  +  ...\ 

JDA  +  (A0  +  Krt  +  tf2r2  +  K3r3  +  A  %  +  ...)  ( - 2 - L - 2_ - 2 - 2 - ] 

+  ^  K0  (R0  +  ATr,  +  K2r2  +  K3r3  +  tf4r4  +  ...)  +  K0  ( R0  +  tfr,  +  tf2r2  +  K3r3  +  K4r4  +  ...) 


</  (R0  +  ATr,  +  A2r2  +  tf3r3  +*%  +  ...) 


j*2/-,2 


+  3  *Va  -  5  6  ■  *3'3+  d  *Vs  "2  d  yi 

*(b)  tt  -j  “» (- ^ *•) =  “ *(Vf *« +  ?) “p (-  <A2) 

Designating r, ,  r2 ,  r3  and  r4  as  drt/dt,  dr2/dt, dr3/dt  and  drjdt  respectively,  collecting  terms  con 
taining  the  same  power  of  K ,  and  rearranging,  we  have 

J  (D*  R0rl  +  JSD  R0)  exp  (-Jf)  dr  =  u(R0^  +  ^  (A3) 


for  AT  and 


Jf  exP  {-st)dt  =  ^ 


r  iri ~  + 


exp  (-sf)  dt 


for  K2  and 


f  R0  r3  exp  (-jf)  dt 

■  jf  }(l - i  {hR‘-  %  '•'fs*)''1 
-  (  -*0  J&Va  *  vi>-(>  -«o  JsVi'1 


iC 


for  K3y  and 

/  R0r4exp(-st)dt  =  f  j-s  ^  rxr2  -  |  ^  rx3  +sr,r2 

*0°  •  -  /7  •  S  3  •  *0*  ,  • 

-  ~D~  W.  +  2  Yd  ~2Dr'  r'-  2D  '<  r> 

+  ^  r»V>  'r3r*  +/?0  ^r2'l-R0  idri2'2 

*  1 §-  *  *•  #  -  r>'>  +  (R°  VJ  - ')  V>| 


exp  (-s/)  dt  (A6) 


for  A4. 

Similarly,  this  process  can  be  continued  to  collect  terms  containing  A5  and  higher.  But  it  can 
be  seen  from  eq  A2-A5  that  the  content,  or  the  complexity  for  determining  r, ,  r2  ...,  increases 
considerably.  As  indicated  in  eq  A6,  in  order  to  evaluate  r4,  terms  of  products  of  rlr2,rtr2,rtr2, 
^tr3'ri,:3*r2^2*ririr2’  ri>':iri'riri  and  rYr2  have  t0  be  determined, and  their  Laplace 
transforms  containing  49  and  1 1 2  terms  with  coefficients  (w  -  1  )3  and  (w  -  1  )4  respectively.  The 
final  Laplace  transforms  of  r.  ,r2,  r,  and  r4  are 
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The  inversions  of  r, , r2, r3  and  r4  are,  respectively, 
r,  =(co-  1) 

r2  =  (co 

,  f2  /8  \  IN  /l3  .8  3 \  Dt 
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Replacing  t  in  terms  of  a  dimensionless  time  r  by  t  •  R\r!D ,  and  substituting  r, ,  r2 ,  r3  and  r, 
into  eq  25,  it  follows 
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Equation  A1 5  can  be  further  simplified  to  contain  only  one  variable  rand  parameters  K  and  to  as 
follows: 

=  A(to  -  1)  (1.1 284  r*4  +  r)  +  A2  (to  -  l)2  (0.7184  r* 

+  0.8634r-  0.3761  r3/2  -0.5  r2)  +  K3  (to  -  l)3 
(0.5555  r*  -0.3451  r  +  1.2628  r3/2  -0.6265  r2-  1.2412  r*/2  +0.5  r3) 

+  K*  [(to  -  l)3  (-0.8165  r*  -  3.2035  r  +  0.2394  r3/2  +  0.0131  r2  +  0.0376  r*/2 
-0.1667  r3  +(to-l)4  (-0.5507  r  ^  +  3.9272  r  +  5.2030  r’^2  -  1 .6403  r2 
-  0.5600  r*/2  +  0.89 10  r3  -  1 .1 660  r7^2  -  0.5  r4)  |  .  (Ah 


For  large  r’s  eq  A16  becomes 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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